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. Abstract Several noncovariant formulations of the electromagnetic self- force 

of extended charged bodies, as have been developed in the context of classi- 
' fH ' , cal models of charged particles, are compared. The mathematical equivalence 

of the various dissimilar self-force expressions is demonstrated explicitly by 
I I deriving these expressions directly from one another. The applicability of the 

■ self-force formulations and their significance in the wider context of classical 

' charged particle models are discussed. 

^ I Keywords self-force • extended charged particle • radiation reaction • 

. classical electrodynamics 
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1 Introduction 



The self-consistent description of charged particles is no doubt the longest- 
standing fundamental problem in classical electrodynamics. The conceptual 
difficulties, which invariably arise from the interaction of the particle with 
its self-generated electromagnetic fields, depend on the particle model that 
is postulated In line with modern particle physics, often a structureless 
• point charge model is adopted, leading to the well-known Lorentz-Abraham- 

Dirac (LAD) equation of motion 0, Q for elementary charged particles. As 
is also well known |^], however, the LAD equation is plagued by irreconcil- 
able deficiencies manifested by runaway solutions, in which the particle mo- 
. mentum grows exponentially toward infinity, or by preacceleration solutions, 

where the particle starts to accelerate even before the onset of any external 
force. A popular remedy is to express in the LAD equation the particle ac- 
celeration perturbatively in terms of the external force Q , which renders the 
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equation stable. However, arguably this procedure does nothing to improve 
the LAD equation intrinsically. It has been put forward that the stabilized 
equation may be considered as more fundamental than the LAD equation 
itself, but the argument usually requires abandoning the point particle limit 
proper by ascribing some structure to the particle Meanwhile, alter- 

native classical electrodynamics of point particles continue to be proposed 



Another well-studied possibility in the development of a consistent charged 
particle electrodynamics is to dispense with the point charge model and 
associated problems altogether, and picture charged particles as extended 
charged bodies. In fact, this was historically the first charged particle model 
that was investigated [14| - |l6| . If the particle size is larger than a very small 
but finite critical length (which is comparable to the classical electron ra- 
dius), the particle dynamics is free of the unphysical runaway solutions fPH - 
[l9l | . This continues to motivate detailed calculations of the electromagnetic 
fields and forces inside accelerated extended charged bodies [l9l - l2^ . More- 
over, extended charged particles keep open the possibility of electromagnetic 
interpretations of inertia (23l |. as originally proposed in pre-relativistic times 
["l5. 16]. Furthermore, the unphysical prediction of singular or preaccelera- 
tion behavior of charged particles at a sudden onset of applied forces can be 
successfully removed by taking into account the finite propagation velocity of 
signals traversing the extended particle [l^l . The latter effect is the classical 
analogue of the concept of " self-dressing" familiar from quantum electrody- 
namics [24] , and plays a role in the measurability of the electromagnetic field 
[25} . Finally, and perhaps most significantly, extended charged particle mod- 
els provide an important way to access point particle models by taking the 
appropriate limit corresponding to vanishing particle radius. Consequently, 
the physical and mathematical consistency of this limiting procedure is itself 
subject of active research [26l - [28j . 



For the above reasons, a thorough understanding of the electrodynamics 
of extended charged bodies is a prerequisite for the development of classical 
charge particle models of both the point charge and the extended charge 
variety. Spread out over a century, however, many formulations of these dy- 
namics have appeared in often dissimilar forms, making a broad comparison 
of results difficult. In this paper, we aim to contribute to a more coherent 
picture of extended particle models by comparing several published expres- 
sions 29-34] for the self- force of a rigid charged body. We demonstrate the 
equivalence of these dissimilar expressions by deriving them directly from one 
another. The self-force F, which is the resultant Lorentz force that is experi- 
enced by the charged body and caused by the self-produced electromagnetic 
fields, reads in noncovariant form 



{pE + Jx B) (fx 



OA 

-P ( V0 + — ) + J X (V X A) 



d^x, (1) 



where p is the charge density, J is the current density, E is the electric 
field, B is the magnetic field, cj) and A are the electromagnetic potentials. 
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Fig. 1 The content of this paper in relation to existing results in literature 



and the integration is over the extent of the charged body. The self-force 
calculations that will be considered here take Eq. ([1]) as a starting point, and 



have the form of series expansions [29 
time [Sll, [S^l , and Fourier integrals 



30.] . definite integrals over retarded 
3J]. Figure [T] shows schematically the 
content of this paper in relation to these publications. It should be stressed 
that the figure represents only a very small fraction of the available literature 
on the subject; correspondingly this paper is not meant as a comprehensive 
review. Rather, the new connections that will be established here complete 
the literature shown in Fig. [U and present the results on a common basis. 
Section [2] gives an overview of the considered existing self-force expressions. 
In section [3] it is shown how the various expressions follow directly from each 
other. Section 0] discusses the applicability of the self- force formulations, and 
puts them in the wider context of classical charge particle models. Like the 
primary equation ([l}, our formulation will be noncovariant throughout. 



2 Existing self-force derivations 

The derivation of the self-force of a rigid charged body requires the evaluation 
of Eq. (ID) by some method. This involves a calculation of the electromagnetic 
potentials appearing in Eq. ([Ij, which necessitates a choice of gauge. The 
potentials f/)*^^^ and A^^^^ in the Lorenz gauge and the vector potential A*^*^^ 
in the Coulomb gauge satisfy the wave equations 

[^^^^')'l'i=^^t)^^'oc^ix,t). (2) 

Here, ^ = cA^^), cA'*^)} and 77 = {cp,J,Jt}, with Jt the diver- 

genceless part of the current density [33| . The scalar potential in the Coulomb 
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gauge is not relevant to the problem, as will be discussed below. In terms of 
the causal Green's function for the wave equation, which equals [s^l 

G(.,.',M') = ^M-^^^, (3) 

\x — x'\ 

with 5 the Dirac delta function, the relevant particular solution of Eq. ^ is 
given by 

<f(a;,i) = ^ jj G{x,x',t,t')n{x',t')d^x'dt'. (4) 

For any given charge distribution and given history of the motion of the 
charged body, n{x',t') is known, so that in principle the potentials can be 
evaluated with Eq. Q, after which the self- force can be determined via Eq. 
(HI). The calculations available in literature where this program is followed 
differ in the order in which the integrations in Eq. (|4]) are carried out. In view 
of the delta function in Eq. Q, it is tempting to start with the integration 
with respect to t' . This immediately yields the well-known retarded integral 
expressions [3^ for the potentials, which indeed are the starting point for 
the self-force calculations presented in sections l2.ll and 12.21 below. However, 
integrating first with respect to t' in Eq. is not the only possibility. For 
certain charge distributions, it is advantageous to start with the integration 
over x', as will be described in section 12.31 Still another possibility is to 
Fourier transform Eq. (|4]), that is, to integrate with respect to the coordinates 
x; this is shown in section [^31 



2.1 Taylor expansion 

Adopting the Lorenz gauge, integration of Eq. ^ with respect to t' yields 
the retarded integral expressions 

Hoc f n{x',tret) ,3 , 



iTT J \X — X'\ 

where now ^ = {(j),cA} and 77 = {cp,J}. In Eq. ([S]), the integration is 
complicated by the fact that 77 must be evaluated at the retarded time 
tret = t~\x ~ x'\ /c, which is different for each volume element cPx'. Jackson 
[34| approaches this problem by expanding 77 in a Taylor series around the 
current time i, 

n=0 ^ ^ 



Substitution of Eq. ^ in Eq. © expresses the potential in terms of quantities 
evaluated at the current time only. Using the result in Eq. (jlj gives, after 
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some manipulations [3J] , the electric part of the self- force 

Here, the integration variables have been changed to r = x — ^{t),r' = 
x' — ^{t) where ^{t) is the trajectory of the center of the charged body, and 
R = r — r' . For a spherically symmetric rigid charge distribution, Eq. ([7]) 
simplifies to 

F = ^ y tlL^^ ff p{r)p{r') \r - r'|""' d'rd'r'. (8) 

In Eq. , the magnetic part corresponding to the last term in square brack- 
ets in Eq. ([IJ has been neglected, so that Eq. ([8]) is the self-force linearized 
in ^ and its time derivatives. In case of harmonic motion ^ = exp(— iwi) = 
^, the series in Eq. ^ can be readily summed, and is proportional to 
exp(za; \r — r'\ /c). Furthermore, writing in Eq. ([5]) the charge distributions 
p in terms of their spatial Fourier transforms, and integrating the resulting 
expression, it is found that [13] 

F = ^ f hm r '"^f .,,, dk. (9) 

Here, the symmetrical convention for Fourier transformed quantities Yj, = 
(27r)^'^/^ J l^(r) exp(— ifc • r)£'r is adopted. The quantity pk is often called 
the form factor of the charge distribution. 



2.2 Lagrange expansion 

As mentioned above, Eq. ([8]) is a linearized approximation to the exact self- 
force due to the neglect of the magnetic term in Eq. ((T]). However, the deriva- 
tion in section [53] is inexact for another reason. Namely, by making use of 
a predefined rigid charge distribution p throughout the derivation (or more 
precisely, using the distribution in the proper frame), it is implied that the 
potentials are generated by a total charge J p{x' ,tret)d'^x' . The latter is in 
general not equal to the true charge of the body f p{x', t)d^x' = g, but rather 
depends on the body's state of motion. To correct for this inconsistency, ei- 
ther the quantity U should be defined in a relativistically covariant way, or 
else the integral ([S]) should be modified to leave the total charge invariant. 
The latter, however, is precisely how the Lienard-Wiechert potentials for a 
moving point charge were devised, as is explained clearly in Ref. d, sec. 19.1]. 
Accordingly, the charged body may be regarded as a collection of infinitesi- 
mal particles moving with the trajectory ^(i) -I- r' and having a fixed charge 
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p{r')d?r' with p the proper frame distribution. The corresponding potentials 
are thus given by 



{c,v} 



poc 

47r ./ R-R-v/c 



(10) 



where v{t) = d^/dt is the velocity of the charged body, and R{t) = x — 
£(t) — r'. An important difference between Eq. ([S]) and the Lienard-Wiechert 
formulation Eq. (jlOp . apart from the different denominator, is that in the 
former the retarded time was known explicitly in terms of the coordinates 
X and x' , while in the latter it is only defined implicitly by the retardation 
condition tret = t — Rjt ret) I c. This complicates the derivation of the self- force 
significantly. Herglotz [23 and Schott [30] proceeded by expanding retarded 
quantities Y in series using Lagrange's reversion theorem l35| . 



Y{tr 



Y{t) 



^ (-1)" 

n—l 



ndY{t) 



dt 



(11) 



Note that differentiation of the quantity i?" in Eq. ((TT|) produces factors of the 
velocity v and derivatives thereof, so that the Taylor series Eq. ([5]) is in fact a 
linearization of Eq. (|TT]) in which all terms nonlinear in v and its derivatives 
have been neglected. Likewise, the potentials Eq. ([5]) are linearizations of the 
Lienard-Wiechert potentials Eq. (flUl) . Now let b be the characteristic size 
of the charged body. On working out the first few terms of Eq. ([TTjl . and 
noting that i? ~ 6 for relevant field points x, it becomes apparent that these 
linearizations are good approximations provided that 



6" d" 



< \v\ 



(12) 



for n > 1 . Roughly speaking, this means that the motion of the body should 
not change significantly on the time scale necessary for light to travel across 
the body, which is the time scale at which self-forces are communicated. This 
condition is known as quasi-stationary motion [s^. It indicates the range of 
validity of the form factor integral Eq. in addition to the condition |t)| ^ c 
associated with the neglect of magnetic forces. 



Substitution of Eq. ([TT|) in Eq. pH)) expresses the potentials in terms of 
quantities evaluated at the current time only. Using the result in Eq. ((1]), and 
performing all integrations, gives a series expansion for the self-force. This 
series has been evaluated explicitly up to cubic terms in the velocity for a 
homogeneously charged sphere of radius b by Herglotz [29] . The linear terms 
are 

6A.og^^ (n + !)(» + 4)(-2Vc)"d"+^| 
7r6 ^ (n-fS)! dt"+2 ' ^ 

n=0 ^ ' 

and dominate the nonlinear terms in case of quasi-stationary motion Eq. 
((T2|l . For a homogeneously charged sphere in rectilinear motion ^(t) = ^(i)ez. 
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Schott [3fl| derived the following closed-form expression including terms up 
to arbitrary order: 



n—O m— 



(m + l)(n + l)(n + 4)(-25)" 
(2m+l)(2m + 3)!(n + 5)! 

a2m+n+2 



(14) 



Qy^2m+n+2 

which reduces to Eq. ([T^ when truncated at m = 



2.3 Direct evaluation 



Sommerfeld [31|, |32] has evaluated the self-generated potentials of a charged 
body by integrating Eq. (U) with respect to x'. For a homogeneously charged 
sphere, this leaves a one-dimensional integral over t' (3^ : 

where Rc{t') = \x — ^{t')\ is the distance to the center of the sphere, and 



X 



4c(t - t')Rc c{t -t') <b-Rc 

9 - [c{t - t') - Rcf b - Rc < c{t - t') < b + Rc . (16) 

c{t -t')>b + Rc 



In virtue of the delta function in Eq. ([4]), times between t' and t' + dt' in 
Eq. (fT5|) correspond to the contribution to the potentials 'F{x,t) generated 
by the charge located within a shell with radius c{t — t') and thickness cdt' 
centered around the field point x. Depending on Rc and t' , this shell may fall 
completely within the charged sphere, or only partially, or not at all, each 
case leading to a different factor x as given by Eq. (ITS)) . Using the potentials 
Eq. Uni) in Eq. ([1]) resuhs in the self-force ^ 




(17) 

in which the integrations are over the time difference t = t ^ t' , and 

G±(r) = [c' - v{t) ■ v[t - r)] tdg{cr±s) ^ dv{t-r)g{cr±s) . 
^ ^ s as s at s 

Here, s = ^(t) — ^(t — r) is the displacement of the charged sphere during the 
time interval r. The upper integration limits in Eq. (jl7p are the roots of the 
equations cr^ ± s(t^) — 2b. These limits demarcate different stages in the 
communication of electromagnetic signals between the parts of the charged 
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sphere that lead to the self-force at the current time t. For subluminal motion, 
the trailing end of the sphere receives electromagnetic signals at time t that 
were emitted by the other parts of the sphere at times between t — r+ and 
t. The signals received by the leading end at time t were emitted by the 
other parts during the slightly longer interval between t — t~ and t. Signals 
emitted at still earlier at times before t — t~ do not arrive at any other part 
of the sphere at time so that the domain r > does not contribute to 
the self-force Eq. ((T71) at aU. 



2.4 Fourier transform 

Bohm and Weinstein (33| have adopted the Coulomb gauge to evaluate Eq. 
(III). The benefit of this gauge choice for the calculation of the self- force of 
a rigid charged body is that the scalar potential (j) equals the electrostatic 
potential corresponding to the instantaneous distribution of charge. Since for 
any pair of charge elements dei and de2 the instantaneous electrostatic force 
on dei due to de2 is the negative of the electrostatic force on de2 due to 
dci, the contribution of (/) to the self- force F integrates to zero identically. 
Therefore only the vector potential has to be taken into account, which is 
given by Eq. Q as before. It can be shown [s^ that a Fourier transformation 
of this equation from the spatial domain x to the wave vector domain fe yields 
the potential 

Ak{k,t) = ^ [ JT.k{k,t)sm[ck{t-t')]dt'. (18) 

Notice that the integration in Eq. ([T5| extends to the upper boundary t, so 
that the potential at time t depends only on currents at past times <' < i, 
that is, Eq. (fT8|) is causal as it should be. Using in the self-force Eq. ([T]) the 
inverse Fourier transform A = (27r)~'^/2 f Ak exp{ik-x)d'^k, and substituting 
Eq. nil), gives [H 

F=—[ dt' [d^k\pk\^exp{ik-s) (19) 

k X [v{t') X k] , v{t) x[kx v{t')] . , 

-T5 COS CKT ; i sin cfcr 

cfc 

Note that the second term in large braces is proportional to and perpendicu- 
lar to the current velocity v{t) of the charged body, and therefore represents 
the magnetic component of the self- force. The first term gives the electric 
component. For a spherically symmetric charge distribution, pk{k) is a func- 
tion of the magnitude of k but not of its direction. In this case, Eq. ([T^ can 
be straightforwardly integrated over angles in fe-space. This reduces Eq. 
to 

j' dt' j^dke\pk\^ J^^(t')_M^^_il£^ jo(fcs)coscfcT (20) 

/ , i\v{t') ■ s\s\ ji{ks) , v(t) X \s X v{t')] , . ■ 

- v{t') ~ ^ ^ ' — ^ ■' ) ^ cos cfcr + —^ ^ —ji (fcs) sin cfcr 

V I ks cs 



eo J- 
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where j„ denotes the spherical Bcssel function of the first kind and order 
n [S^l- In Eq. (^0]) . the first two terms in large square brackets represent 
the electric component of the force and are given in Ref . fssf ; the last term 
gives the magnetic component. The integral over k containing Bessel function 
kernels has the typical form of an inverse Fourier transform in spherical co- 
ordinates [38|. In section [21 we will derive the other self- force representations 
given in sections 12.11 to 12.31 from this Fourier integral. 



3 Equivalence of the self-force expressions 

3.1 Fourier integral and integral over time 

Sommerfeld derived for the piecewise function x given by Eq. (|16p the integral 
representation [32| 

86^ f°° h{kb) smjkR,) sm[ck{t - t')] 

^ = V X k 

Substituting this representation in Eq. p5p . and using the result in Eq. ([IJ, 
gives the self-force [32] 



9q^ 



dT I dk[ji{kb)f (22) 



I g s 
— — \v(t — T)jo(ks) sin cfcr] — — v(t) ■ v(t — r)] — ji(fcs) sin cfcr 
k ot ^ ■' s 

Performing the integration over k indeed yields the force Eq. ([T7)) . We now 
show that Eq. (|22l) is equivalent to Eq. (1^0)) that was derived by Fourier 
analysis of the potentials in the Coulomb gauge. Note that the integrands 
of both equations already have a similar structure due to the form of the 
integral representation Eq. ([?!]) . Performing the differentiation d/dt in Eq. 
([^ using the property 

gives, after changing the integration variable back to t' and rearranging, 
.2 r* 



9q^ 



27r2eofo2c J_ 



dt' I dk [ji{kb)f sinck{t - t') 







—niks)^ (23) 



1 dvit') f {[vit)~v{t')]-s}vit') [vit)-v{t')]s \ . 

Next we integrate by parts the first two terms in the large square brack- 
ets with respect to i', choosing for the differentiated factors respectively 
fi(t') = sji(ks)/s and f2(t') — jo{ks) sinckr. To carry out this integration 
unambiguously, it is necessary to replace the lower integration limit t' = —oo 
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by t' — —a, and take the limit a ^ oo afterwards. With the help of the 
relations 

a/i fShiks) . \ k[v{t')-s]s vit') . 



dt' \ ks J s 

= — ji(ks) sin cfcr — ckjo(ks) cos cfcr, 

or s 



the resulting self-force is 

0„2 



27r2eo&2 

where 



poo 

lim / [n{kbf{B + I)dk, (24) 



B = 



I = 



v(t') . s 1* 

— ; — jo (ks) sin cfcr — - — ji (ks) cos cfcr 

— ^ ^ '^^ — L j jp (ks) cosckr 

a L ' 



v{t') - i — — )^ — cos ckT+^ ^ —ji {ks) sin ckr 



dt'. 



Taking in Eq. (j24p the limit a ^ oo oi I presents no difficulties, and yields 
precisely Eq. (^0]) . specialized to a homogeneous sphere that has the form 
factor 

Therefore Eq. (IM)) is equivalent to Eq. (pn| . provided that the boundary 
term B vanishes. This can be shown to be the case as follows. B evaluated at 
t' = t vanishes since sinc/cr = and s{t') — at t' — t. In the limit t' — oo, 
the first term of B is zero trivially when v{~oo) — 0. When v{—oo) ^ 0, 
it must be that s{t') oo and hence jo{ks) as t' ^ —oo, so that the 
first term does not contribute in this case either. The second term of B, on 
the other hand, vanishes at t — 5- — oo only when v{—oo) ^ 0. Namely, when 
v{—oo) = it is possible that s(— oo) = S has a finite value. In that case, 
the boundary term makes a contribution to Eq. (1241) proportional to 

/"OO 

lim / h{k) cos ck{t ~ a) dk, (26) 

where h{k) = Ijijkb)]^ jiikS)/k. However, Eq. ([25)) evaluates to zero by the 
Riemann-Lebesgue lemma jS^. Hence = for all possible v{—oo), so that 
the force Eq. ([24)) is indeed identical to the force Eq. (l20l) that was derived 
by Fourier analysis of the potentials in the Coulomb gauge. 
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3.2 Fourier integral and form factor integral 

As discussed above, the self-force Eq. ^ in terms of a form factor integral 
is valid for quasi-stationary motion Eq. (IT^ and \v\ <§; c, and for the spe- 
cial case of harmonic motion. In order to compare Eq. Q with the self-force 
derived in section [2^ the latter should be specialized accordingly. This may 
be effected by expanding the integrand of Eq. ([T^ in a Taylor series around 
t' = t, and linearizing the result by neglecting all terms nonlinear in v and 
its derivatives. The extremely involved full expansion, in which all nonlinear 
terms have been kept, is given in Ref. [2l[. Formally, such use of a Taylor 
series to describe the integrand on the infinite interval — oo < t' < i is ques- 
tionable because the series may have a finite radius of convergence. However, 
for subrelativistic motion electromagnetic signals are communicated between 
parts of the charged body on a time scale ~ b/c, so that only the small in- 
terval t~b/c<t'<t significantly contributes to the integral in Eq. ([T^. 
This can be seen by noting in Eq. (1191) that the integrand only contributes in 
the domain |fc| < because the form factor \pk\ ~ elsewhere, and that 
the integral over this domain averages out due to the sinusoidal functions 
unless ckr < 7r/2, that is, unless t — b/c ^t'<t. Proceeding on this basis by 
Taylor-expanding, neglecting nonlinear terms, and integrating over angles in 
fe-space, yields 



Writing T2"coscfcT = {-!) " c-^^'{d/dk)^" cos ckr and 

t2"+i cos cfcr = (-l)"c-2»-i(d/dfc)2"+i sin cfcr, as is suggested in Ref. [2ll |. 
and integrating by parts with respect to k repeatedly, gives 




(27) 




(28) 




in which 




2n-l 




2n 
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All boundary terms Bn'^ and Bn'^ vanish identically. At fc = cx), this is 
because pfe(oo) = for any finite charge distribution. At A: = 0, the terms 
with odd m are zero because \pk\ is an even function, and those with even 
m vanish because they contain sincA;T as a factor. The quantity ^^j'Kln^ = 
j'^'' can be interpreted as the symmetric cosine transform of the function 
j2n(fc) = {dj dk)'^'"'k^ \pk\', likewise, ^j2pKln^ = Jn^ is the symmetric sine 
transform of j2n+i{k)- Therefore the double integrals in Eq. (pS)) reduce to 
the single integrals 

(30) 




j2n+lik) 



Here, the identity k ^ sin kxdk — 7r/2 and the variable x = ct have been 
used. With the help of Eqs. (EH) -([301), the force Eq. ^ reduces to 



3e^^ c2"(2n)! 

n— 



(2n+l)c V7o A: / di^n+i 

This expression now has the manageable form of a series in terms of the 
derivatives of the current velocity, with coefficients that are readily calculated 
from the form factor of the charge distribution. In the next section, we will 
specialize this result to a homogeneously charged sphere, and show that it 
is equivalent to the series expansion Eq. (fT3|) obtained by application of 
Lagrange's reversion theorem. Here, we apply Eq. (1311) to the case of harmonic 
motion, for which v — —iuj^o exp{—iujt) = —ioj^. Since {d/dt)^v = {—iuj)^v, 
Eq. ((3T|) then becomes the sum of two ordinary power series in the quantity 
Lj I c. The series corresponding to the first line of Eq. (1311) may be interpreted 
as the even part of the Taylor series of the function p(k) = around 
K = 0, evaluated at k — lo/c. Similarly, the series in the second line may 
be identified with the odd part of the Taylor series of p{k) around k — k, 
evaluated at k = k + uj/c. Summing these two series therefore results in 

F=^j(^IJpik + u^/c)-p{k-cv/c)]f+'-^[pic./c)+p{-u^/c)]y (32) 

Noting that p{k) is an even function, the integral in Eq. ((32|) may be rec- 
ognized as the Hilbert transform of p{k) in a less common notation (39l |. 
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Accordingly, by changing variables it may be shown (39l | that 

f . !l^ff/"|lW^£,,+ „ (33, 

where j denotes the Cauchy principal value. Here, the second line of Eq. ([5^ 
has been interpreted as a residue. Eq. ([55]) is identical to the force Eq. © 
derived from a Taylor expansion of the retarded integrals for the potentials. 



3.3 Fourier integral and Lagrange expansion 

In the previous section, we derived the series expansion Eq. (|3T]) that ex- 
presses the linearized self-force in terms of the derivatives of the current 
velocity of the charged body, for a general spherically symmetric charge dis- 
tribution. We will now specialize this result to a homogeneously charged 
sphere, and show that this yields the self-force Eq. ((T3)) that was obtained 
from series expansion of the retarded potentials. Evaluation of Eq. ([5T|) using 
the form factor of a homogeneous sphere Eq. (|25p requires determination of 
the quantities 



x=0 
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The first of these equals (2n) ! times the coefficient of x^" in the Taylor series 
of [ji (x)]'^ around a; = 0. By squaring the ascending power series of the Bessel 
function [13|, it is thus found that 



" ^^m!(n-m- l)!(2m + 3)!!(27i-2m + l)!!' ^ ' 



Writing factorials in terms of Pochhammer symbols {p)q = T{p + q)/r{q) 
with r the Gamma function [s^l, Eq. ([55]) becomes 

'"-r(|)r(n)r(n + i)ir'o ' ^''^ 

Here, it has been used that {p)-q = (— 1)'^/(1 — p), The series in Eq. 
defines a Gauss hypergeometric function with unit argument [S^] ■ Evaluating 
this hypergeometric function, and converting Gamma functions to factorials, 
results in 

(n + l)(n + 2)(2n + l)- ' 
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Establishing T„ is more involved. The squared Bessel function [jiix)]'^ can 
be expanded in a series of Bessel functions with doubled argument [41[ . This 
gives 



2m + 2 



m=0 



(2m+ l)(2m + 3) 



(38) 



where J denotes the cylindrical Bessel function of the first kind (33]. The 
factor X in the denominator can be removed with the help of the recurrence 
relation 2pJp{z)/z — Jp_i(z) + Jp+i(z). Subseq uently, the integrand of r„ 
in Eq. p4p is found by application of the expansion [42| 



dPJgjz) 

dzP 



u=0 



Jq-p+2u{z)- 



This gives 



oo 2n4-l 

EE 



(2m + 2)(-l)'' 
(2to+ l)(2m + 3) 



2n + 1 

u 



(39) 

Co^4m-2n+2ii(2a;) 

(40) 



m— li— 

+ C'2^4m-2n+2«+2(2a;) + C4 J4m_2n+2M+4 (2a;) + CqJ im.-2n+2u+e{'^x) 

where 

Co 



1 



C2 
C4 



(4m + 3)(4m - 2n + 2u + 1) ' 
1 

(4m + 5)(4m - 2n + 2m + 5) ' 
2(4m + 4) 



(4m + 3) (4m, + 5) (4m - 2n + 2-u + 3) ' 
2(4m + 4) 



(4m + 3)(4m + 5)(4m - 2n + 2w + 3) ' 

Substituting this expansion in Eq. ([M]). the integral T„ can be evaluated 
trivially because Jq°° Jp{z)dz = 1 for arbitrary p > —1 [37| . Therefore T„ is 
given by Eq. if each Bessel function is replaced by 1/2. The remaining 
double series can be summed in closed form. The sums over u of the various 
terms have been tabulated [4^; together they evaluate to 



7r(-l)"(2n + 1)! 
= 4r(n+i)r(n+|) 



(41) 



E 



(1) 



1 

m V 4 



-!(! + t)™(^ + t) 



The series in the second line of Eq. (j41l) is derivable from the series 



Un{z) - E 



(1) 



1 

m V4 



m— 



-!(l + t)™(¥ + t). 



(42) 
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which defines the generahzed hypergeometric function [42 

" 1 1^2k 3~2n 

^ A A 



Un{z) - ^F^ 

Comparing Eqs. (|1T|) and (li^ . it is found that 



ll+27t 9+27t 
4 4 



7r(-l)"(2n + 1)! 



Arin 



\)r{n. 



AUn{l)+ Un{z)D{z)dz 



with D{z) = - z^/^. The integral in Eq. (gH) is equal to [i] 



1 1 l-2rt 3-2k 
-^2 4 4 

5 11+271 9+2n 
2 4 4 



(43) 



(44) 



(45) 



Eqs. (I43p - (|45p define T„ in terms of two generalized hypergeometric functions 
of unit argument; for both functions closed form expressions in terms of 
Gamma functions exist [i^ . Writing these expressions in terms of factorials 
yields, after considerable reduction, 



7r(2n + l)(-4)" 
(n + l)(2n + 3)(2n + 5)' 



(46) 



Finally, having the quantities Sn and T„ at our disposal, the self-force Eq. 
dSD) can be evaluated. Combining Eqs. and dM]) gives 



neocb'^ ^ 



(2n)(2n + 3)(-26) 



2n 



1 (ff-'^V 



(2n + 4)! c^" c;i2" 
(2n+ l)(2n + 4)(-26)2"+i 1 



(47) 



J2n+1. 



(2n + 5)! 



r.2n+l ^^2n+l 



Taking the two terms in large braces together by relabeling the summation 
index gives precisely Eq. (fT3)) . It has thus been shown that the self- force 
obtained by Fourier analysis of the potentials in the Coulomb gauge is equiv- 
alent to the force derived by Lagrange expansion of the potentials in the 
Lorenz gauge. 



4 Discussion and conclusions 

In this paper, we have demonstrated the equivalence of a number of pub- 
lished expressions for the self-force of a rigid charged body. These included 
the form factor integral Eq. ([9]) of Jackson [s^l based on a Taylor expansion 
technique, the series Eqs. and of Herglotz and Schott [3fl] re- 
sulting from a Lagrange expansion, the integral over retarded time Eq. (jl7p 
of Sommerfeld [33] obtained by direct integration of the fundamental equa- 
tion, and the Fourier integral Eq. ((20|) of Bohm and Weinstein [s^ derived 
by Fourier analysis adopting the Coulomb gauge. To this list we may add 
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our result Eq. (j3ip . which is intermediate between an integral and a series 
representation. 

The various expressions differ in their degree of generality. The Fourier 
type results of Jackson and of Bohm and Weinstein allow calculation of the 
self-force for arbitrary spherically symmetric charge distributions, whereas 
the direct integration and Lagrange expansion methods, which stay in the 
spatial domain throughout, require specialization to a particular charge dis- 
tribution at an early stage. In addition, several approximations implied in 
the derivations may be discerned. The neglect of magnetic forces in Eq. ([7]), 
and of nonlinear terms in the series Eq. (jl3p , restrict the validity of the corre- 
sponding self-forces to nonrelativistic velocities. More subtle is the improper 
use of the charge distribution at the retarded time, as discussed at the start 
of section 12.21 which further constrains the results based on Taylor rather 
than Lagrange expansions to the quasi-stationary motion defined by Eq. 
(|12p . This restriction does not apply to the more correct self- force expansion 
Eq. (jl4p . On the other hand, numerical evaluation of this series by trun- 
cation is only possible when the series converges sufficiently rapidly, which 
is when (6/c)("+iH^^/'^*)^""^^^C < (Vc)"(d/di)"C, a condition even stronger 
than quasi-stationary motion. In contrast, the Fourier technique leading to 
Eq. (1^ does not involve the explicit evaluation of retarded source densities 
in the first place, so that the issue is avoided altogether. Therefore, Eq. (^0]) 
may be regarded as the most general self-force expression in this paper, in 
the sense that this result does not neglect any magnetic or other nonlin- 
ear terms, and is not restricted to quasi-stationary motion or any particular 
charge distribution. At the same time, Eq. (j20p is suitable for numerically 
stable evaluation as it does not involve truncated series. 

As emphasized in the introduction of this paper, the electrodynamics of 
extended charged bodies is important in relation to point charge models of 
elementary particles, since the two concepts are connected by an appropriate 
limiting procedure. It is therefore illustrative to see how some well-known 
point particle results follow from the self-force expressions discussed in this 
paper. This connection is most straightforward in case of the series expansion 
Eq. p3p . which is valid for nonrelativistic velocities. In the point particle 
limit 6 — > 0, the first terms become dominant and evaluate to the well-known 
nonrelativistic self-force [4, chap. 21] 

AUes dv (Pv 

with Ues = 3(?^/(207reo6) the electrostatic energy. The first term on the right 
of Eq. (|48l) has the appearance of an inertial term; correspondingly the quan- 
tity 4L/es/(3c^) is usually interpreted as an electromagnetic contribution to 
the mass of the particle. However, the factor 4/3 violates relativistic mass- 
energy equivalence, and its removal by more sophisticated arguments has 
been the subject of much discussion [l3, [H, |46| - |48| . The inertial term is 
commonly absorbed in an effective mass of the particle, leaving the second 
term on the right of Eq. (|48| as the radiation reaction force usually featuring 
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in the equation of motion of the classical point particle, in the nonrelativistic 
limit. The derivation based on Taylor expansion in section I^TTI produces the 
limit Eq. (^5]) as well. The first two terms of Eq. ^ are 

3c2 \87re JJ |r-r'| J dt' ^ ' 

1 / . / . , 9 , \ d V 



The quantity in large parentheses in the n = term is just Ues, so that Eq. 
(HH]) equals the inertial term of Eq. (|^5)) . The quantity in parentheses in the 
n = 1 term is the total charge, and Eq. (|50| reduces to the radiation reaction 
term of Eq. (j48p . Moreover, these observations are independent of the partic- 
ular charge distribution one chooses to model the particle with 0, chap. 21]. 
Less obvious is the point particle limit of the Fourier integral representation 
of the self- force Eq. (1^0]) . which, as we have shown, is equivalent to the series 
Eq. (fT3)) . and should therefore reduce to Eq. (|48)) as well. However, this is 
difhcult to see from the Fourier integral directly. In contrast, from our inter- 
mediate result Ec^. (pij) the limit Eq. is readily obtained. Namely, the 
term n = in the first line of Eq. pil) vanishes for any finite pk- The term 
n = in the second line gives, after integration by parts, the force 



2 dv f \pk\'^ 



Seoc^ dt J 



d^k. (51) 



By writing out pk as the Fourier transform of the charge distribution p(r), 
and simplifying the result by using the identity j exp \—ik ■ (r — r')] k^^d'^k — 
2Tr^/ \r - r'l [Si, Eq. is reduced to Eq. and hence to the inertial 
term of Eq. (|48p. The radiation reaction force, in turn, is generated by the 
term n = 1 in the first line of Eq. (1311) . Exploiting in the latter spherical 
symmetry by writing pk — \/2/t: p{r)ji){kr)r'^dr^ performing the differ- 
entiations with respect to k and evaluating the result for fc = 0, yields the 
force Eq. ((5(7)) . None of the other terms of Eq. ([HT]) contribute in the point 
particle limit. This is because pk will tend to a constant as the support of 
p[r) shrinks to a point, and therefore the differentiations with respect to k 
will make these terms vanish. This confirms that also the Fourier integral 
Eq. (PT|) correctly reduces to the well-known self- force Eq. (H51) . 



An important aspect of the rclativistic dynamics of an extended body, left 
out of the discussion so far, is the way is which the body maintains or changes 
its shape while being accelerated. On the one hand relativity theory requires 
a velocity dependent Lorentz contraction; on the other hand instantaneous 
contractions are impossible due to the finite velocity at which information 
about velocity changes can propagate through the extended body. Acceler- 
ated motion, therefore, necessitates some notion of relativistic rigidity, the 
precise formulation of which is still being studied [10, [Mi . Often the approx- 
imation of Born rigidity [s^ is applied, valid for adiabatic velocity changes, 
where it is assumed that the accelerated body always maintains its shape in 
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its continuously changing instantaneous rest frame [53|. The corresponding 
Born rigid extended electron model was advocated by Lorentz [l^ . A signif- 
icant strength of this model is that it provides a natural way to cure the 4/3 
problem described above, by assuming a negative pressure inside the elec- 
tron that balances the Coulomb repulsion of the distributed charge [i^ . The 
work done by this pressure during Lorentz contractions removes the factor 
4/3 precisely [l^[23- (Another possibility is to adopt a manifestly covari- 
ant definition of four-momentum of systems involving electromagnetic fields 
[iJjSZe^Sj)- However, none of the self- force derivations discussed in t his p a- 
per corresponds to the Lorentz model. Rather, the Abraham model is 
implied, in which the charged body maintains its shape in the frame of an 
observer at rest. Differences between the self- forces according to each model 
do not yet appear in the linear terms, however, so that the limit Eq. (|48)) 
is model independent. In the Lagrange series Eq. (fH|) . on the other hand, 
the sum over m of the nonlinear terms n — indeed agrees [30| with the 
model of Abraham and not with that of Lorentz. Similarly, the time integral 
Eq. (fT5)) . and therefore also the equivalent Fourier integral Eq. ((20|) . coincide 
with the Abraham model [3l|. Since the contracting Lorentz model seems 
more in line with special relativity, it would be valuable for classical charged 
particle theories to derive a Born rigid version of the Fourier representation 
Eq. ((20|l of the self-force. A covariant formulation based on formal cut-off 
procedures has been developed in this direction Yet, it should be real- 
ized that it has not been established that a consistent extended body model 
of elementary particles, if any, should necessarily respect rigidity. Therefore 
alternative postulates, such as the Abraham model or other non-rigid models 
[131, may continue to prove their value. What is more, the current state of 
technology is starting to enable experimental conditions in which the elec- 
tromagnetic self-force of macroscopic charged systems, such as high-density 
electron bunches [s^ and ultracold plasma bunches [s^, becomes significant. 
It would be interesting to see to what extent the self-force formulations in 
this paper can model these evidently non-rigid systems. 
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